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INTRODUCTION 


In  a  recent  paper  (ref  1)  we  presented  quantitative  results  of 
diffraction  model  (ref  2)  calculations  of  the  temperature  dependence  of  the 
resistivity  of  amorphous  magnesium  zinc  (a-MgZn).  The  phase  shift  expansion 
of  the  scattering  matrix  elements  as  derived  by  Evans  et  al  (ref  3)  was  used 
with  the  phase  shifts  adjusted  to  give  the  observed  magnitude  of  the 
resistivity  and  to  satisfy  the  Frledel  sum  rule.  Saturation  effects  (ref  4) 
were  included  by  invoking  the  Pippard-Ziman  mean  free  path  constraint  (refs 
1,2,5)  on  the  electron-phonon  interaction.  A  single  effective  scattering 
potential  was  used  in  place  of  the  full  alloy  scattering  matrix  elements. 
Comparison  of  these  computed  results  with  the  data  of  Matsuda  and  Mlzutani 
(ref  6)  in  a-MgZn  showed  remarkable  agreement  between  theory  and  experiment  on 
all  the  many  features  of  the  temperature  dependence  of  the  resistivity.  It 
was  therefore  considered  worthwhile  to  refine  the  calculations  by  using  the 
alloy  scattering  matrix  elements  with  phase  shifts  computed  for  Mg  and  Zn  in 
an  extended  version  (ref  7)  of  the  Ziman-Faber  substitutional  model  (ref  8). 
Selected  results  pertaining  to  a-MgZn  are  presented  here.  Additional  results 
relating  to  low  resistivity  alloys  in  general  are  given  in  Reference  9. 


THEORETICAL  MODELS 


We  employ  a  form  (ref  7)  of  the  Ziman-Faber  liquid  alloy  substitutional 
model  (ref  8)  which  is  appropriate  for  an  amorphous  Debye  solid  (ref  10).  The 
result  for  the  electrical  resistivity  is 


12  An  IKK 
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References  are  listed  at  the  end  of  this  report. 
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where 


|U(K)|2  -  c^InCK)  -  t2(K)|2lP(K) 

+  |c!t!(K)  +  c2t2(K) | 2SP(K) 

and  and  t^(K)  are  the  concentration  and  t-matrix  for  the  ic^  component. 


The  resistivity  structure  factor  is  given  by 


SP(K)  -  e"2w(K)a(K)  + 


3R  K 


4  .A  .2 


4  d“q 


+  -  ( - ) 2  /  d(--)(— )  n(x)  (n(x)+l)F(qA)  f  ---  a(K+q) 

T  2kF  0  qD  qD  0  4  w 


where  a  ■  3(2hkF)2/Mkg9^. 


The  corresponding  cross  term  integral  is  given  as 


,  K  .2  4  ..q  _q  .2 


IP(K)  -  e-2wCR)  + - (-—)  f  d(--)(-~)  n(x)[n(x)+l]F(qA) 

T  2kF  0  qp  qp 


The  partial  geometric  structure  factors  in  the  substitutional  models  are 


assumed  to  be  identical  and  are  given  by 


1  ■*>  ♦  ♦ 
a(K)  "  -  l  exp[iK.(m-n) 
N  m,n 


where  m  and  n  represent  the  average  ion  positions. 


The  scattering  matrix  elements  for  the  ith  component  is  given  by 


2irh3  r  .  ihfcJ(EF) 

tj(K)  - - - - I  (2£+l)8in  nfcJ(EF)  P^cos  9) 

m(2mEF)l/2n0 


where  n£j(EF)  are  the  scattering  phase  shifts. 


The  Pippard-Ziman  constraint  (ref  5)  on  the  electron-phonon  Interaction 


reflects  the  breakdown  of  the  adiabatic  approximation  due  to  finite  electron 


mean  free  paths  A  and  is  incorporated  here  through  the  Pippard  function  (ref 


11)  F(qA)  where 
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(7) 
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2  y  tan“*  y  3 

F(y)  -  -  [ - t - 1 

w  y-tan-1  y  y 


In  Che  long  mean  free  path  limit  F(qA)  ♦  1  and  Che  above  expressions  reduce  Co 


standard  diffraction  model  resulcs. 


The  remaining  Cerms  in  Che  above  expressions  have  their  usual  meaning; 


ocher  details  and  definitions  are  given  in  Reference  9. 


RESULTS  AND  DISCUSSION 


The  phase  shifts  for  these  calculations  were  obtained  by  constructing 


muffin  tin  potentials  with  the  Mattheiss  prescription  (ref  12)  using  Herman 


and  Sklllman  (ref  13)  neutral  atom  wave  functions.  For  Mg  the  Kmetko  (ref  14) 


value  of  0.75  of  full  Slater  exchange  was  used.  For  Zn  the  value  of  0.85  of 


full  Slater  exchange  was  used  since  this  places  the  d-levels  at  approximately 


9  eV  below  Ep,  which  is  appropriate  for  Zn  alloys.  The  phase  shift  values 


were  then  computed  using  modified  versions  of  computer  routines  given  by 


Loucks  (ref  15).  The  i  ■  0  to  t  ■  3  values  at  Ep  in  the  72  percent  Mg  alloy 


for  Mg  are  -0.175,  0.085,  0.034,  and  0.001,  respectively.  For  Zn  the 


corresponding  values  are  0.354,  0.294,  -0.057,  and  0.002. 


The  other  parameters  used  in  the  computations  were  based  on  the 


following:  x-ray  data  (ref  16)  for  a-Mg7Zn3  show  the  main  peak  in  the 


structure  factor  Is  located  at  kp  *  2. 7 A-1  for  the  a-MgZn  composition  range. 
We  assume  that  kp  varies  so  that  2kp/kp  ■  1.11  is  appropriate  for  the 


composition  range;  the  average  ionic  mass  was  used  for  M;  a  Debye  phonon 


spectrum  Is  assumed  and,  since  z  ■  2,  we  took  qp  -  kp;  specific  heat  data 


(ref  17)  give  6q  ■  295K  for  the  low  temperature  limit,  but  we  selected 


SB 


iw; 


200K  for  the  resistivity  Debye  temperature  to  fit  the  room  temperature  TCR. 
(For  crystalline  Mg,  0D  ■  390K  and  Sr  -  340K,  while  for  Zn,  Sq  *  310K  and  Sr  ■ 
175K  (ref  18).)  A  geometric  structure  factor  of  Percus-Yevick  form  with 
packing  fraction  of  0.525  is  appropriate  for  amorphous  alloys  (ref  18)  and  was 
used  here. 

Detailed  experimental  data  for  a-Mg]_xZnx  for  .225  <  x  <  .325  were 
obtained  by  Matsuda  and  Mizutani  (ref  6).  All  compositions  show  negative 
TCR’s  at  the  high  temperatures  and  resistivity  maxima  ranging  between  40  to 
60K;  most  also  show  broad  shallow  minima  centered  at  approximately  10K. 

Another  feature  found  by  Matsuda  and  Mizutani  is  that  the  resistivity  varies 
as  (T-Tm)^/^  for  x  >  T^»  where  Tj(  is  the  temperature  at  the  resistivity 
maximum . 

The  results  of  the  present  calculation  are  shown  in  Figures  1  and  2.  The 
computed  magnitude  of  p  is  about  40  uftcm.  In  Figure  1  we  see  the  p  vs.  T 
behavior  computed  from  the  standard  diffraction  model  (Eqs.  (1)  through  (7) 
with  F(qA)  -1):  a  maximum  and  a  negative  TCR  at  room  temperature  are 
obtained  in  the  calculation,  however,  TM  is  higher  than  observed  by  about  a 
factor  of  two,  the  calculated  magnitude  of  the  TCR  at  room  temperature  is  too 
small  by  at  least  a  factor  of  four,  and  the  size  of  the  maximum  is  too  large 
by  two  orders  of  magnitude.  The  low  temperature  T2  limit  occurring  in  the 
standard  diffraction  model,  (i.e.,  p  -  p0(l+BT2)  with  B  s  SxlCT6^2)  holds  up 
to  about  15K.  This  is  not  consistent  with  the  data  either  since  the  mean 
value  for  B  obtained  by  Matsuda  and  Mizutani  (ref  6)  is  0.5  x  L0_6K”2  and  the 
temperature  range  for  this  approximate  quadratic  dependence  is  considerably 


extended. 


Resistivity  vs.  temperature  curves  for  32. 5%  Zn  composition 
dp  •  p(T)  -  p(300);  no  saturation. 


Figure  2  shows  the  result  with  saturation  effects  included.  The  value  of 
qDA  was  set  at  11.6  to  correspond  to  the  observed  p  value  at  x  ■  .325  and  to  a 
saturation  resistivity  of  200  pflcm.  The  previous  calculations  (ref  1)  with 
the  effective  potential  approximation,  including  saturation,  reproduced  many 
of  the  observed  features  such  as  the  presence  of  minima,  the  correct  magnitude 
and  approximate  average  position  of  the  maxima,  the  (T-T[j)'^  dependence,  and 
the  reduced  magnitude  of  B.  The  present  results  show  these  features:  for 
example,  the  value  of  B  for  x  ■  .325  is  approximately  0.5  x  10-6K-2  with  an 
extended  T2  range.  In  addition,  the  observed  shape  and  position  of  the  minima 
are  obtained;  the  trends  in  size  and  position  of  the  minimum  and  maximum  in 
resistivity  and  the  magnitude  of  the  resistivity  are  also  consistent  with  tf 
data. 

CONCLUSIONS 

The  results  underscore  the  point  that  effects  of  finite  electron  mean 
free  paths  are  important  in  the  temperature  dependence  of  the  resistivity  of 
amorphous  alloys,  even  in  alloys  with  resistivities  as  low  as  50  u£5cm.  Only 
qualitative  agreement  with  experiment  is  obtained  if  saturation  is  ignored. 
Earlier  work  (ref  2)  showed  that  not  even  qualitative  agreement  with  the 
diffraction  model  could  be  obtained  at  higher  resistivities  unless  these 
effects  are  considered.  The  detailed  agreement  between  the  present 
computations,  which  include  saturation,  and  the  many  features  of  the  p  vs.  T 
data  of  a-MgZn  is  remarkable. 
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ATTN:  TECH  LIB  -  DRXMR-PL 
WATERTOWN,  MA  OI272 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 
ATTN:  CHIEF,  IPO 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  2 


DIRECTOR 

US  NAVAL  RESEARCH  LAB 
2  ATTN:  DIR,  MECH  DIV 

CODE  26-27,  (DOC  LIB) 
WASHINGTON,  D.C.  20375 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
1  ATTN:  AFATL/DU 

AFATL/DUG 
EGLIN  AFB,  FL  32542 


COMMANDER 

US  ARMY  HARRY  DIAMOND  LAB 
ATTN:  TECH  LIB 
2800  POWDER  MILL  ROAD 
ADELPHIA,  MD  20783 


METALS  &  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  LAB 
505  KING  AVENUE 
COLUMBUS,  OH  43201 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CTR 
ATTN:  TECHNICAL  LIBRARY  1 

CODE  X212 

DAHLGREN,  VA  22448 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  SMC AR-LCB-TL , 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


NO.  OF 
COPIES 


1 

1 


1 

1 


1 


